URC97064 

Turbulent Distortion of Condensate Accretion 

R. Hazoume J. ChabiOrou,** J. A. Johnson III 
CeNNAs, Florida A&M University, Tallahassee, FL 32310 (USA) 

Abstract 

When a simple model for the relationship between the density-temperature fluctuation correlation and mean 
values is used, we determine that the rate of change of turbulent intensity can influence directly the accretion rate of 
droplets. 


Considerable interest exists in the accretion rate for condensates in nonequilibrium flow with icing’ and the 
potential role which reactant accretion can play in nonequilibrium exothermic reactant processes. 2 Turbulence is 
thought to play an important role in such flows. It has already been experimentally determined that turbulence in- 
fluences the sizes of droplets in the heterogeneous nucleation of supersaturated vapors. 3 This paper addresses the 
issue of the possible influence of turbulence on the accretion rate of droplets. 


According to a nonequilibrium model developed earlier 4 the droplet growth rate for flow with changing local 
pressure is given by: 


dr(t ) 
dt 
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dt 



(i) 


with B= (2s)/[r c T(R/m n )] wheres = surface tension coefficient r c = condensate density; and m n = molecular weight 
of vapor. In this approach, we use the hardcore model for a liquid droplet in which dN, the number of simple mole- 
cules of water in the shell between r and dr, is given by dN =(k/r”)(4rcr 2 )dr. N n is a correction coefficient such that 

0<N n ^4 and N n =[(n-3) n -%(n-4) n - 2 ] for n>4. P^ is obtained through the Clausius-Clapeyron equation. For adia- 


batic flow of a compressible ideal gas the mean value of the square of the local speed of sound is defined as 
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b 2 = 3 = =y— ( 2 ) 

( d P) 7T*_ P 

where: g = Cp/Cv is the ratio of the specific heats; g is the mean value of the condensate mass fraction of water; P 


is the mean pressure function; S is the mean entropy; and p is the mean density function. 


We can readily write from (2) that 

P P 

— r = cnst and — ~ = cnst (3) 

p (P)(r=i) 

Equations (2) and (3) were written assuming there are fluctuations of the principle physical properties around their 
mean values as, e.g., g = g + g’ with g' = O. However g' 2 #0. This means that we can consider 

P — P + P' , p — p + p' , g = g + g' , and T — T + T' . The equation of state 

p = p(i— ^ + ^lZ1) RT (4 ) 

Pi Pv 

becomes _ 

P + P' = (p + p')(}—^P- + (0 ?~ 8 ~ 8 )R(f + T) (5) 

Pi P v 

with 0) o =the initial humidity, p i =the molecular mass of the carrier gas, and R=the ideal gas constant. Averaging 
equation (4) yields: _ 

p = (Ll^o + ^o_8 )R( fp- + If 7 ) (6) 

Pi Pv 

where a new function, p ‘T’, has appeared, this is the density-temperature correlation function involving the fluc- 
tuations of density and temperature. We can use this result in eq. (5) above. 
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Now we make the important simple assumption 

W = -*PT. ( 7 ) 

in which eq. (7) is a definition of X; with this definition, X is without consequences for any specific experiment. The 
equation of state then becomes 

4 = /(g)(l-A)*f = — (8) 

P 7 


where 


m= 1 ^ 


(o Q -g 


(9) 
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We must also make it clear that the physical properties of interest are functions of space and time. That is, 

T = T(x i ,t),P =P(x i ,t),p = p(x i ,t),g =g(x i ,t),X = X(x i ,t) (lo) 


where i = 1,2,3, xi is the i component of the space variable, and t is the time. To use equation (10), we can use as the 
simplest adequate physically reasonable model for A(xi,t) to be — A(jt { - , t ) = CXp— ft. The exponential char- 
acter of the decay is consistent with the usual long-term behavior of turbulent correlations. With this model, the 
density temperature correlation function depends only on time and obeys a decay law which is exponential. Finally, 
it is then usual (with U and L as velocity and length parameters respectively) to express all these relations in a non- 
dimensional form: vi =Ki/U; yf=x(/L; and t=tU/L. 


We can now include the influence of fluctuations in temperature (as well as the implicit fluctuations in P and 
p) in a turbulence distorted refinement of the local mean square of the local speed of sound a 2 (y 2 , t) by explicitly 

P(y , t) a 2 

requiring _ 2 ’ — = — . This form is reconciled with equation (8) above by 

P (y 2 ,i) 7 
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We have shown in equation (3) that Pa(b 2 ) <k7 ~^ . Thus 
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in which/is the real positive constant introduced in the definition of 1. We finally obtain from equation (1) 


dr(r) 6(7 ( n - 3)" 3 y / _dLoga 2 (y„ r) 

dT p c (—)T 4 ) n " 2 (7 - 1) exp/r + 1 dT 


( 13 ) 


When we define the instantaneous turbulent intensity as proportional to d 2 {y 2 ,r), then the interpretation of 

equation (13) is straightforward. (i)At small values of t and constant turbulent intensity, ( da 2 {y 21 r) / dr = 0 ), the 
rate of change in droplet size is determined entirely by the density correlation exponential factor/ . (ii)At large val- 
ues oft and constant turbulent intensity, ( da 2 (>’ 2 , t)/ dr =0), the rate of change in droplet size has no dependence 
on turbulence at all. (iii)If the turbulent intensity is not constant, i.e., ( da 2 (y 2 , T )/ dr # 0 ), then the rate of change 

in droplet size is sensitive to the turbulence at all times. The case da 2 ( y 2 , f)/ dr ^0 requires that a measure of the 
actual evolution in the strength of turbulent fluctuations be determined and that a model for determining the connec- 
tion between the mean square local speed of sound and the evolution of turbulence be determined. 
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The overall features just described are summarized in Fig. 1. It has already been shown^ that, at fried ther- 
modynamic conditions, the size of a droplet in condensing flow is influenced by the strength of the local turbulence 
and, by inference, the value of the Reynolds. Equation (13) above sets conditions for the sensitivity of the conden- 
sate accretion to turbulence, based on basic physics principles, which can now be tested by experiments. 

This research was supported in print by NASA Grants NAGW-2930 and NAG 2-291 and by DOE Contract 
DE-FG04-95AL87309 to Florida A&M University. 
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Figure 1. Sensitivity in the Evolution of Turbulent Distortion of Droplet Growth Rate with the Persistence of the Turbulence Pulse. 
The symbol/is defined in the text; T 2 is the laboratorytime; T, is the persistence time of adistortion in the turbulence; and XT is the 
product of the left hand side of eq. (13) with the prefactors of the right hand side's bracketed expression. The plots show the 
differences between a short puke, in the top display, and a long and relatively unchanging turbulence in the lower display. 
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